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Three BCS superconductors Sa, Sb, and S and two short normal regions Na and in a three- 
terminal SaNaSNfoSfo setup provide a source of nonlocal quartets spatially separated as two correlated 
pairs m Da and Sfe, if the distance between the interfaces NaS and SN?, is comparable to the coherence 
length in S. Low-temperature dc transport of nonlocal quartets from S to Sa and can occur in 
equilibrium, and also if Sa and S^ are biased at opposite voltages. At higher temperatures, thermal 
excitations result in correlated current fluctuations which depend on the superconducting phases (j)a 
and (/)b in Sa and Sb. Phase-sensitive entanglement is obtained at zero temperature if Na and 
are replaced by discrete levels. 

PACS numbers: 74.50.+r,74.78.Na,74.45.+c,03.67.Bg 



Regarding the manipulation of entangled states, quan- 
tum nanoelectronics is on the way to address the same 
fundamental issues with electrons as quantum optics does 
with photons. An entangled quantum state has a den- 
sity matrix distinct from that of any "hidden-variable" 
theory. Two-particle entanglement can be probed [1] via 
the violation of the Bell inequality [2]. Multiparticle en- 
tanglement also has a high potential, for instance it Q 
can be used to implement error correction codes. 

Concerning superconductivity, two-particle entangle- 
ment can be generated at normal metal-superconductor 
NSN interfaces, by extracting a split Cooper pair from 
the BCS condensate of electron pairs [4-7]. We show in 
this Letter that a nanoscale three-terminal superconduct- 
ing setup can produce nonlocal quartets separated as two 
pairs in different electrodes, therefore opening a route for 
a new generation of entanglers which could be controlled 
by an electromagnetic field. One must stress that here 
quartets are absent in the bulk superconductors which 
instead carry ordinary BCS pairing. This is in contrast 
with the destruction of the "ordinary" pair condensate 
in certain arrays of Josephson junctions [8|, and its rela- 
tionship with topological quantum computation 0, ^| . 
Microscopically, nonlocal quartet transmission appears 
here as a generalization of so-called Cooper pair splitting 
at a double NSN interface, and it can be characterized by 
interference and noise. Recall that an Andreev pair in a 
normal metal electrode N results from the emission of a 
charge 2e from S at a NS interface, by the process of An- 
dreev reflection (AR). At a double NSN interface, spin- 
entangled [11] or energy-entangled [l^ pairs can be pro- 
duced through crossed (or nonlocal) Andreev reflection 
CAR [iMsi involving evanescent quasiparticle states in 
S, on the coherence length ^s- CAR coexists with normal 
transmission through S without electron-hole conversion 
(elastic cotunneling EC) [16]. CAR or EC can be se- 



lected by their different Coulomb interaction energy [I3 , 
by their spin sensitivity [gI, fl6j, by their distinguishing 
energy dependence [18], or by their different signature in 
the nonlocal conductance and in the zero-frequency shot 
noise cross-correlations [H, [ll|, [l9|, I1q| . The new effects 
considered here do not require more advanced technology 
than the experiments on split pairs already realized with 
metallic structures [13] or with quantum dots [14]. 

In this Letter, a route to the production of nonlocal 
quartets is proposed on the basis of bunching of two An- 
dreev pairs in a superconducting beam splitter made of 
conventional BCS superconductors. Indeed, two Joseph- 
son junctions separated by a distance ds of the order 
of the coherence length ^5 of S can be coupled by non- 
local coherent effects [21]. Here we study microscopi- 
cally all the possible nonlocal effects and discuss their 
physical consequences. Remarkably, in a three- terminal 
SaNaSN5S5 structurc, nonlocal quartets can be sepa- 
rately transmitted as two correlated pairs in Sa and Sfo. 
Non-local quartet transmission proceeds through double 
crossed Andreev reflection (dCAR), which coexists with 
double elastic cotunneling dEC. The latter process yields 
Cooper pair transmission between Sa and S5 [22]. The 
phase-sensitive dCAR is a new coherent nonlocal quan- 
tum channel, which has no direct analog for incoher- 
ent multiple Andreev reflections \2^. The elementary 
charges involved in dCAR and dEC are doubled as com- 
pared to CAR and EC. The four processes of CAR, EC, 
dCAR and dEC (see Fig. [Tj) will be treated on an equal 
footing [24], as well as AR at each of the SaNaS or SN5S5 
interfaces which transfers Cooper pairs between S and Sa 
or S5. 

Cooper pair splitting in a SaNaSN^S^ structure with 
arbitrary interface transparency can be described by gen- 
eralizin g th e Andreev-Kulik-Saint- James bound states 
(ABSs) [25|42S|- Those states are coherent superpositions 
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FIG. 1: (Color online). Panel (a) shows a SaSSb structure, 
where Na and have not been represented for clarity. The 
processes taking place in three-terminal SaNaSN^Sb [panel 
(b)] and Sa^SpS^ [panel (c)] structures are as follows: dou- 
ble crossed Andreev reflection (dCAR, red long dashed lines) 
producing a nonlocal quartet (a spatially separated pair of 
pairs), double elastic cotunneling (dEC, dotted black lines) 
exchanging pairs between Sa and S^, and crossed Andreev re- 
flection (CAR), thermally activated above the gaps of Sa and 
Sfe. In addition, elastic cotunneling (EC) and local Andreev 
reflection (AR), not shown in the figure, also take place in 
SaNaSNfoSfe. With suitable gate voltages, the n- and p-doped 
semiconductors on panel (b) have a vanishingly small density 
of states at negative and positive energies, respectively. They 
filter CAR and dCAR, and eliminate EC, dEC, and AR. 



of electrons and holes, forming in a short single channel 
junction a doublet at opposite energies {E-^E^\ and 
carrying Josephson currents in opposite directions. In 
the SaNflSN^Sfo setup, the phases of Sa, S5 and S being 
(\)a^ (\>h and (\)s^ and taking the bias voltages Va, V5 to be 
zero, the total currents J^, Ih are given at zero temper- 
ature T = by la[<i>aAh) - (2e/^)En=12(^^n/S(/>a), 

h{<i>aAh) - (2e/^)En=i,2(^^n/S(/>6). The are the 
energies of the ABSs formed by hybridizing the ABSs 
of both junctions through dCAR and dEC. They de- 
pend on both phase differences (5^^ = (\>a — (\>s and 
= — (j)s. If ds ^ £,s^ this leads (at lowest or- 
der in the dCAR and dEC processes) to la = I^i^^Pa) + 
I^^^{5(^a - Hh) + /^^^^((^(/>a + Hh) and h = itiHh) - 
I'^^^{S(l)a-S(l)b)^I'^^^^{S(j)a^S(j)b). Production of nonlo- 
cal quartets (dCAR) and pair transmission (dEC) couple 
the coherent dc Josephson currents in Sa and S5 by the 
inverse crossed inductances {L~^)a,b — S^a(^a, ^5)/S^5 
and (L~^\^a = dh{4^ai(t^h) / d(j)a^ which is an extension 
of the concept of crossed conductances in a NaSN^ struc- 
ture. 

Voltages Va^ V5 are applied now on the electrodes 
Sa, S5 {Vs = is the reference voltage). Yet a dc 
Josephson current can flow from S to Sa and S5 if 
ds [29!], in addition to the standard ac Joseph- 

son currents. The intensity of this current can be seen 
as a synchronization of the phases of the ac oscillations. 
Indeed, considering for simplicity low transparency con- 



tacts, the double Josephson junction is described by the 
Hamiltonian H = Hioc + HdCAR + HdEC + 2e{hsVs + 
fiaVa^hbVb), with H\oc = -Ej[cos{S(l)a{t))^cos{S(l)b{t))], 
and HdCAR = -Ef^^cos{SMt) + ^bW), HdEC = 
—Ej^^cos{S(l)a{t) — 6(j)i}{t)). The currents are obtained 
from Hamilton equations, for instance: 



T^dCAR 



Ia{t)=Ejsm[5Mt)]- 



(1) 



^^^^^ sm[SMt) + Hb{t)] + sm[Ha{t) - SMt)]' 

Applying opposite voltages Va = —Vb = V leads to a dc 
Josephson effect for nonlocal quartets because the phase 
combination 5(j)a{t)-\-S(l)b{t) = e{Va-^Vb)t/h-\-S(j)a-\-S(j)i) is 
then time independent. The Josephson effect for nonlocal 
quartets becomes ac only if the energy eVa-\-eVb acquired 
by the quartet when separately transmitted into Sa and 
Sb is finite. This result holds for any transparency. 

The dc Josephson effect for nonlocal quartets is fur- 
ther considered for a SafiSpSb junction biased at oppo- 
site voltages, where the previous Na and N5 metals have 
been replaced by n- and p-doped semiconductors. The 
conduction band edge on one side (n type) and the va- 
lence band edge on the other {p type) are at zero energy. 
The gaps in the density of states of the n- and p-doped 
semiconductors filter the processes with positive energies 
in n, and with negative energies in p [l^] [Fig.[T]^c)]. This 
excludes both the local Josephson effect and the nonlo- 
cal dEC, thus leaving at T = only the nonlocal dCAR 
as a coherent coupling between the condensates [13] • In 
this ideal situation, one obtains a perfect superconduct- 
ing beam splitter operating at the scale of the coher- 
ence length, and producing correlated pairs of Cooper 
pairs flowing in the leads Sa and Sfo. Notice again that 
the biases Va and should be opposite in the coherent 
Josephson regime, while they are equal in a normal beam 
splitter NaNNfo or a Cooper pair splitter NaSN^, where 
quasiparticles are emitted instead of pairs. 

Let us now discuss the noise cross-correlations in a 
SaNaSN^Sfo structure. Zero- frequency thermal noise is 
present in the absence of applied voltage for sufficiently 
transparent single junctions [3l|,l32|. Finite values are ob- 
tained for all the components of the correlators Sij{t) = 
{SIi{t + t')5Ij{t))^ where 5Ik is the current fluctuation 
in Sfe {k = a, 6). This equilibrium noise, due to ther- 
mally activated fluctuations between ABSs carrying op- 
posite currents, is phase sensitive because the popula- 
tion of quasiparticles in thermal equilibrium exchanges 
charge with the condensate. The thermal noise can be 



very large for a long inelastic lifetime of the ABSs [31, 
CAR and EC lead to components S^^^ and 5'^^ of Sa,h, 
which are independent of 5 (pa and ^^5. They correspond 
to CAR and EC assisted by thermal activation over the 
gap A (Fig. [Tj). On the contrary, dCAR and dEC re- 
sult in thermal fluctuations between hybridized ABSs. 
dCAR corresponds to random emission and absorption 
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of nonlocal quartets between S and S^, S5. dEC corre- 
sponds to random transmission of pairs between and 
S5. The contributions 6'^^^^ and S'^ec of dCAR and 
dEC to Sa,b depend, respectively on the phase c ombina - 

Generalizing Refs. 



tions 5(j)a-\-^4^b and 8(j)a — S(j)i). 
the zero-frequency noise cross-correlations Sa,h in the ab- 
sence of applied voltage is written as 



1 



— V — 



(2) 



where the lifetime of the Andreev states l/rj shows the 
relevance of pair-breaking effects in the superconduc- 



tor [2l| . Such noise correlations could be large at temper- 
atures T* comparable to the ABS energy level difference. 
The crossover temperature T* is strongly reduced as the 
interface transparency increases. 

Calculations of the equilibrium nonlocal inverse 
inductance (L~^)a,5(^^a, ^^b) and cross-correlations 
5*0, 6 (^^a, ^^^5) are based on microscopic Nambu-Keldysh 
Green's functions [29| in which the interfacial hopping 
amplitude is accounted for by a self-energy. Arbitrary 
values of the temperature and of the normal-state trans- 
mission coefficient Tat can be treated in equilibrium be- 
cause the time convolutions in the Dyson equations then 
simpHfy into products of Green's function depending only 
on energy. The microscopic calculations carried out for 
a three-dimensional ballistic superconductor apply to a 
voltage range in which the proximity effect is negligi- 
ble. Sa^h goes to zero at T = because it is ther- 
mally activated over the ABS gap. {L~^)a^h saturates 
at low T to the zero-temperature response of the con- 
densate. As seen from perturbation theory in the tun- 
nel amplitudes, both dCAR and dEC contribute with 
a positive value to Sa,b{^4^a = 0,J^5 = 0). Thus, 
for tunnel contacts, Sa,b{^(l>a = 0,(5^5 = 0) is positive 
at the smah T/A = 0.1 [see Fig. Efb)]. The oscilla- 
tions in —{L~^)a,b{^<i)a^^4^b) [Fig. EJa)] match those of 
Sa,b{^<i)a^^<i)b) (scc Fig.[2)3). They reflect the distinguish- 
ing phase dependences of dCAR and dEC. Increasing Tn 
has the effect of favoring the transmission of pairs by 
dEC from Sa to S5 (or from S5 to Sa), and disfavoring 
their transmission as a pair of holelike quasiparticles by 
dCAR. The cross-shaped variations of S<pb) 
and Sa,b{^4^a^^4^b) for intermediate Tn can be under- 
stood from the bound state energy level minima in the 
{S(l)a,S(l)b) plane for 8(j)a,Hb ^ {0,7r,27r}. 

This superconducting beam splitter also generates en- 
tanglement between pair numbers in the two branches, 
as pairs are emitted two by two. Indeed a split Josephson 
current due to nonlocal quartets connects coherently pair 
number states \Na)\Ns)\Nb) with states \Na + 
2)|A/'5 + 1), entangling the added pairs in Sa and 5'^ [29[. 
To illustrate this, let us replace the N junctions by quan- 
tum dots Da, D5. Indeed phase-sensitive entanglement is 
obtained in a SaD^SD^S^ structure biased at voltages 
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FIG. 2: (Color online). The plots show — (L ^)a,b [panels 
(a), (c), (e)] and Sa,b [panels (b), (d), (f)] as a function of 
{6(l)a,S(l)b) for T/A = 0.1, Tiv = 4.10"^ [panels (a),(b)], Tn = 
0.64 (panels c, d) and Tn = 1 [panels (e), (f)]. The ratio 
A/ As = 0.1 is used. 



Va = -Vb = V larger than A, with As > A, eV. 
In a first step, D^ and D5 are supposed to carry spin- 
degenerate orbitals. In addition, the gate voltages are 
such that Da (D5) has a level at eVa {eVb) but no level 
at —eVa {—eVb). Quasiparticles are transmitted from Sa 
to Da (from Sb to D5) but local Andreev processes be- 
tween Da and S (D5 and S) are not possible. In the 
limit of large gaps, three terms contribute to the effec- 
tive Hamiltonian of the two coherently coupled levels at 
energies eVa (in Da) and eVb (in D5): (i) The exchange 
of pairs between Sa and Da (Sb and Db): y^AR,a{b) = 
-a[exp(i0a(6))c^(5) .^c+^)^^ + /i.e.]. (ii) Cooper pair split- 



ting: TicAR 



+ ^+ . 



h.c.]. (iii) Produc- 



tion of nonlocal quartets: HdCAR 

h.c] . Exact diagonalizations of HAR,a + ^AR,b + ^car 
^dCAR result in an entangled ground state characterized 
by a positive concurrence, which depends on the values 
of 5(j)a and 5(j)b via the combination 5(j)a + ^(j^b typical 
of dCAR [2^. The Coulomb interaction Hamiltonian is 
Uu = Uha(b),t'^aib),ii with ha(b),a the number of spin-a 
electrons in dot a (in dot b). As U increases, the zero- 
temperature concurrence of the ground state remains fi- 
nite because of virtual excitations coupHng to dCAR. 

Entanglement can also be more directly assessed by 
showing that no classical correlation can account for 
the crossed noise of nonlocal quartets. Let us consider 
SaDaSD5S5 or SapSnS5 setups, biased at opposite Va and 
Vb larger than A, and with A^ A, eV [assumption 
(Al)]. The notation {Na{t,r))p stands for the average 
number of electrons transmitted into electrode Sa in the 
time interval [t^t + r]. The average over the hidden- 
variable density matrix p = J dXpa{X) PbW is noted 
(• • • )p [11]. We make the assumption (A2) that each 
subsystem a and b is separately described by quantum 
solid-state physics: only communication through S cor- 
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responds to a hidden variable. The assumption A2 leads 

to {SNait,T)SNbit,T))gu = {SNait,T)SNbit,T))p, with 

{dNa{t,T)5N,{t,T))p (3) 
= J d\f{\){5Na{t,T))^{5Nb{t,T))^, 

where /(A) is the probabihty density of the hidden vari- 
able A. The value of SNa{t^r) for the specific value A of 
the hidden variable is denoted by {6Na{t,T))x' 

An additional assumption (A3) concerning the setup 
is made in order to simplify the discussion: the gap 
is larger than the bandwidth W of the superconductors 
Sa and S5, and the linewidth broadening in S is 7^5 = 0. 
The three assumptions (Al), (A2), and (A3) imply that 
charge transport is blocked at any temperature because 
the Da and D5, or the n and p energy filters [assumption 
(Al)], suppress all Josephson processes taking place lo- 
cally within each subsystem a or b for any realization of 
the hidden variable A [assumption (A2)]. No quasipar- 
ticle is transmitted from S to Sa (from S to S5) within 
subsystem a (b) for any value of A (assumption A3). The 
equality {5Na{t^r)5Ni){t^r))qu = is then obtained be- 
cause {6Na{t,r))x = {dNb{t,r))x = 0. The following 
two statements are in conflict: (i) The cross-correlations 
are finite; and (ii) The nonlocal processes is a classical 
communication (related to some hidden variables) rather 
than the quantum mechanical CAR, dCAR. 

For the Sa^iSpS^ structure considered above, the cross- 
correlations are vanishingly small at zero temperature 
and the Bell-like argument does not imply entanglement 
in this case. Indeed, the ground state of a Josephson 
junction is entangled only if the total number of pairs 
is fixed However, the denomination phase- sensitive 
entanglement is appropriate at zero temperature for the 
SaDaSDfoSfo structure because the cross-correlations are 
finite at T = 0, and the Bell-like argument is in agree- 
ment with the direct calculation of the concurrence, being 
also finite. 

To conclude, we propose to generate nonlocal quar- 
tets in the solid state, by producing correlated pairs in a 
superconducting beam splitter involving three supercon- 
ductors. The nonlocal inductance and the phase-sensitive 
thermal cross-correlations may be probed in future exper- 
iments. Of particular interest is the possibility to obtain 
entanglement at zero temperature from current-current 
cross-correlations in a SaD^SD^Sb structure. Non-local 
Shapiro steplike experiments are also promising for in- 
vestigating dCAR and dEC, in the sense that resonances 
due to dCAR or dEC could be obtained if the ac voltage 
oscillations are synchronized in two incoming channels Sa 
and S5. 
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